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Abstract 

Resonant P-wave interactions can be described by a minimal zero-range model defined by a 
truncated effective range expansion, so that the only 2-body interaction parameters are the inverse 
scattering volume 1/ap and the P-wavc effective range rp. This minimal model can be formulated 
as a local quantum field theory with a P-wave interaction between atom fields and a molecular 
field. In the two-atom sector, the model is renormalizable, but it has unphysical behavior at high 
energies, because there are negative-probability states with momentum scale rp. In the sector with 
three atoms, two of which are identical, renormalization in some parity and angular-momentum 
channels involves an ultraviolet limit cycle, indicating asymptotic discrete scale invariance. The 

— 1/3 

Efimov effect occurs in the unitary limit ap , rp — >■ 0, but this limit is unphysical because there are 
low-energy states with negative probability. The minimal model can be of physical relevance only 
at energies small compared to the energy scale set by rp, where the effects of negative-probability 
states are suppressed. 

PACS numbers: 03.65.Ge, 21.45-v, 34.50.Cx 
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I. INTRODUCTION 



Systems with an S-wave resonance near a threshold arise in many fields of physics, including 
atomic, condensed matter, nuclear, and high-energy physics. They are important, because 
they exhibit universal behavior that depends on the energy of the resonance, but is otherwise 
insensitive to the nature of the constituents and the form of the interactions between them, 
as long as they have short range. The universal aspects of the few-body problem have been 
studied in detail pj. The universal properties depend on the S-wave scattering length a. If 
the Efimov effect occurs in the three-body problem, the universal properties also depend on a 
three-body scale parameter upon which observables can only depend log-periodically. In the 
many-body problem with identical fermions with two spin states, the universal properties 
have also been studied thoroughly, theoretically and also experimentally using ultracold 
atoms [2H5] . If the Efimov effect occurs in the three-body problem, the many-body problem 
is more complicated and may not even be well-defined. 

Systems with a P-wave resonance near a threshold also arise in many fields of physics. 
A prominent example in nuclear physics is the low-energy P-wave resonance in n— ^He scat- 
tering which corresponds to the unstable ^He nucleus [6]. The universal aspects of such 
systems are not as clear-cut as in the S-wave case. They involve at least two parameters: 
the P-wave scattering volume ap and the P-wave effective range rp. The many-body physics 
of fermions with two spin states that interact through a P-wave resonance has been studied 
using mean-field approximations [7HTT]. If the P-wave interactions are strongly resonant, the 
many-body physics may reveal qualitatively new features that are not captured by mean-field 
approximations p!2] . 

Ultracold atoms provide a promising laboratory for the experimental study of both few- 
body and many-body systems of particles with a P-wave resonance near threshold. The 
inverse scattering length ap can be controlled experimentally and can be made arbitrarily 
large by tuning the magnetic field to a P-wave Feshbach resonance. The first such studies 
were for a P-wave Feshbach resonance between identical fermions using atoms [T3] . 
Fermionic ®Li atoms have been studied near a Feshbach resonance between the lowest two 
hyperfine spin states [T5] . The binding energies of P-wave dimers and the atom-dimer 
and dimer-dimer inelastic collision rates have been measured IT7] . P-wave Feshbach 
resonances have also been observed in a mixture of fermions and ^^Rb bosons [IH] . One 
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obstacle to studying the dependence on ap is that P-wave Feshbach resonances in ultracold 
atoms are usually very narrow. 

In this paper, we study the renormalization in the three-body sector for a minimal field- 
theoretic model of a system with a P-wave threshold resonance. In Section |TT| we explain 
why the model can only be physically relevant at energies small compared to the energy 



scale set by rp. In Section III, we write down the integral equation that can be used to 



calculate the bound-state spectrum of triatomic molecules in the model. In Section |IV 
we identify the angular-momentum and parity channels of the three-atom sector in which 
a three-body parameter is required by renormalization. Finally, in Section |V], we discuss 
whether the renormalization properties of the minimal model can be relevant to any real 
physical systems. 

II. TWO-BODY PROBLEM WITH A P-WAVE RESONANCE 

We begin by discussing P-wave scattering at low energy E = A;^/(2/i), where /i is the 
reduced mass. If the particles interact through a short-range potential, the P-wave phase 
shift has an effective range expansion in powers of k"^ that has the form 

cot Sp{k) = -1/ap + lrpP + ... , (1) 

where ap is the scattering volume and rp is the P-wave effective range, which have dimensions 
(length)^ and l/(length), respectively. If the only important interactions are in the P-wave 
channel and if the effective range expansion is truncated after the second term, the scattering 
amplitude reduces to 

/c^ COS 

fk{0) = -—. rp; ■ (2) 

— 1/ap + -^rpk'^ — ik-^ 

An equivalent expression is obtained if the interactions are dominated by a coupling to a 
P-wave Feshbach resonance: 

cos 

^"^^^ ^ iiy-ky2fi)/g^-tk^ • 

The parameter rp = which is negative definite, controls the strength of the coupling 

to the Feshbach resonance. The detuning of the resonance from the threshold is controlled by 
the combination 1/aprp = fiu. We will consider P-wave interactions for which the scattering 
amplitude has the minimal form in Eqs. (|2| or ([s]) with the two parameters ap and rp. 



A point in the parameter space that is of particular interest is ap = rp = 0, which 
is called the unitary limit because the unitarity bound for scattering is saturated at this 
point. Since they provide no length scale, the interactions are scale invariant in the unitary 
limit. It is easy to see that the unitary limit cannot be realized with interactions through 
a short-range potential. If the interaction potential vanishes outside the range R, there is 
an upper bound on rp called the Wigner bound [I9]. For |ap|^/^ ^ R, the Wigner bound 
reduces to rp < —2/R. If we try to take the zero-range limit i? — )■ to justify the truncation 
of the effective range expansion in Eq. ([T]), rp is driven to — oo. The constraint rp < — 2/i? 
can also be derived by demanding that the probability for a bound state to be in the region 
r>Ris less than 1 [201 EH- 

The behavior in the zero-range limit is improved if the potential with range R is sup- 
plemented by a van der Waals tail that falls off like —Cg/r^. In this case, the expansion 
in Eq. ([T]) also includes a linear term bpk [22]. In the zero- range limit i? — ?■ 0, the coeffi- 
cients bp and rp are determined by the scattering volume ap and the van der Waals length 
Pq = {2fiCe/h'^y^^. If we also take the limit \ap\^^^ > /^g, bp ^ and rp -3A/(3q. Thus 
the low-energy expansion of fkid) reduces to Eq. ([2]), but the scale of rp is set by 1//3q. The 
truncation of the expansion is justified only at low energies \E\ ^ 

Real atoms interact through a short-range potential with a van der Waals tail, but they 
can also have couplings to diatomic molecules. The scattering volume ap can be controlled 
and made arbitrarily large by tuning the magnetic field to a P-wave Feshbach resonance 
where one of the molecules crosses the two-atom threshold. The conditions for fk{0) to be 
well-approximated by the effective range approximation in Eq. ^ have been studied by 
Zhang, Naidon, and Ueda [23] . 

Nishida has pointed out that the assumption that the scattering amplitude has the simple 
form in Eq. ([2]) or ([s]) up to arbitrarily large momentum k necessarily implies the existence 
of states with negative probability [21]. In the case of a short-range potential, the problem 
is related to the Wigner bound. If the bound rp < —2/R is violated, the probability for 
a bound state to be in the region r > R exceeds 1 and that excess probability must be 
cancelled by negative probability from the region r < R [2^ . 

The problem of negative-probability states can also be seen directly from the expression 
for the scattering amplitude in Eq. It can be expressed as a function of the energy 
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E = e/2n: 



fk{9) 



2iJ,E cos 6 



(4) 



l/ap + ^irpE - (-2/iE - ief/'^ ' 



The poles in E of this function are the energies of bound states. We set l/ctp = for 
simphcity. If rp < 0, as required by the Wigner bound, the denominator then has zeroes at 
-E = and E = — rp/8yU, which correspond to one bound state at threshold and another with 
binding energy rp/8/i. The scattering amplitude can be expressed as a sum of contributions 
from the two poles and a function of E that is regular at the poles: 



Note that the residues of the two poles have different signs. The two bound states will con- 
tribute to unitarity sums with the opposite sign of the residues in Eq. ([s]). The probabilities 
of the bound states with binding energies and rp/8/i are positive and negative, respec- 
tively. If rp > 0, there is only a negative probability pole with energy on the physical 
sheet. A negative-probability state with energy is fatal for any physical interpretation of 
the threshold region. If rp < 0, the negative-probability state with binding energy rp/S/x 
is not necessarily fatal, because its effects are suppressed at sufficiently low energy. Thus a 
sensible physical interpretation requires rp < and \E\ ^ rp/8/x. If l/ap is nonzero, the 
analysis is more complicated. We refer to a pole as shallow if it has energy \E\ < rj^/Sfi and 
as deep if it has energy \E\> Vp/Sfj,. The different cases can be classified as: 

1. l/ap < 0: There is only one deep pole on the physical sheet. It has negative proba- 
bihty. 

2. < l/ap < |rp|^/54: There are two poles on the physical sheet, one shallow pole with 
positive probability and one deep pole with negative probability. 

3. l/ap > |rpp/54: There are two poles with complex energies and complex residues 
on the physical sheet. This case violates standard analyticity assumptions for the 
S-matrix and should be discarded |25] . 

The problem of atoms whose pair interactions give the scattering amplitude in Eq. ^ 
can be formulated as a quantum field theory (for a review of field theoretical models of 
atom-atom interactions see e.g. Ref. [26]) with two scalar fields ipi and 1IJ2 (which annihilate 
atoms of types 1 and 2, respectively) and a vector field d (which annihilates a diatomic 
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molecule). We take the masses of the atoms to be mi and m2. The only interaction is 
a P-wave contact interaction that allows transitions between the diatomic molecule and a 
pair of atoms of types 1 and 2. The form of the Lagrangian is constrained by Galilean 
invariance. It consists of kinetic terms, an energy offset for the diatomic molecule, and a 
P-wave interaction term: 

^-pi (4 + ^ '"^ ■ (4 " w^>^") " 

+Aod^ -d-gofi [dt • (V^2^VV^i/mi - (iVz^a)^^!^) + (h.c.)] . (6) 



The parameters Aq and Qq are bare parameters that depend on the ultraviolet cutoff. The 
parameter t] can be chosen by normalization of the field d to be either +1 or —1. If 
?7 = +1, d is an ordinary field which, in the absence of interactions, would annihilate an 
ordinary positive-probability molecule. If 77 = —1, d is a ghost field which, in the absence 
of interactions, would annihilate a negative-probability molecule. 

The two-body problem for atoms 1 and 2 can be solved analytically. The minimal model 
defined by the Lagrangian in Eq. (|6| is renormalizable in the two-atom sector [6J . The bare 
parameters go and Aq can be tuned as functions of the ultraviolet momentum cutoff A in 
such a way that the scattering amplitude reduces to Eq. ^ in the limit A — > 00. The 
renormalization conditions can be written 

— = — Ao + — A^ 7a 
ap g^n Svr 



rp Sn 2 



1^ - (7b) 



In order to have finite limits as A — )■ 00, the two terms on the right sides of both Eqs. (7a) 



and (7b) must have opposite signs. A finite limit for ap requires that Aq be negative and 
that Ao/gQ scale like A^. A finite limit for rp requires that gQ scale like A~^ and that 
?7 = — 1. Thus d must be a ghost field. If there were no interactions, it would annihilate a 
negative-probability molecule. To see whether there are any negative-probability particles 
in the presence of interactions, one needs to examine the scattering amplitude in Eq. ([2]). As 
we have seen, rp < is necessary to avoid negative-probability particles near the two-atom 
threshold. The energy restriction \E\ <^ rp/8fi is then necessary to ensure that unphysical 
effects associated with negative-probability particles are suppressed. The simplest observable 
that can go wrong is the two-body bound-state spectrum. It contains unphysical low-energy 
states if the energy restriction is not observed. 



Nishida also presented an algebraic argument for the existence of states with negative 
probability in this model [21]. In the unitary limit Op = rp = 0, the model has scaling 
symmetry as well as Galilean symmetry. It therefore also has nonrelativistic conformal 
symmetry [27]. The asymptotic behavior of the scattering amplitude in Eq. (|2| at large 
k implies that the dimer field d has scaling dimension 1. However Nishida and Son have 
used the nonrelativistic conformal algebra together with the assumption that all states have 
positive norm to prove that primary operators can only have scaling dimensions greater than 
or equal to 3/2. The violation of the bound in this model implies that there must be states 
with negative norm. 

III. P-WAVE STM EQUATION 

We now consider the three-body problem with resonant P-wave interactions between pairs 
of atoms. For three identical bosons, there can be no P-wave interactions. The case of three 
identical fermions has been studied thoroughly by Jona-Lasinio, Pricoupenko, and Castin 
|21j . They considered a two-channel model with atoms and a diatomic molecule that have 
interactions with a finite range h. They found that rp has to be negative and that there is a 
lower bound on \rp\ that is proportional to 1/6. In the three-fermion sector, they calculated 
the spectrum of triatomic molecules (trimers), atom-dimer scattering lengths, and three- 
body recombination rates. For \ap\ ^ 6^, the trimer spectrum at energies small compared 
to h'^/mb'^ can consist of either or 1 trimer with positive parity and either or 1 trimer 
with negative parity. Three-body recombination was also studied previously by Suno, Esry, 
and Greene using the adiabatic hyperspherical representation of the Schrodinger equation 
for three identical fermions interacting thorough a short-range potential that is tuned to 
give a large P-wave scattering volume |28| . 




In the next simplest case of a three-body problem with resonant P-wave interactions, 
there are two identical atoms of type 2 and a third atom of type 1. Atoms of types 1 and 
2 interact through a P-wave resonance and we assume that the interaction between the 
identical atoms can be neglected. We describe this system using the model defined by the 
Lagrangian in Eq. ([6]). Transition amplitudes for three atoms can be calculated exactly nu- 
merically by solving a single-variable integral equation that is analogous to the Skornyakov- 
ter-Martyrosian (STM) equation for S-wave contact interactions [29]. The equations for 
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different orbital angular momentum and parity quantum numbers can be decoupled. 
The bound-state problem reduces to solving an eigenvalue equation for the energy E < Q 
that is a homogeneous integral equation of the form 



B{p) 



3f-l 



iJ+i 



V- 



dqR' {p,q,E)D{q,E)B{q), 



(8) 



7r(2J + 1) 

where the prefactor z/ is +1 or —1 if the identical atoms are bosons or fermions, respectively. 
The upper limit A of the integral is an ultraviolet momentum cutoff that should be taken 
to oo. The dimer propagator D{q, E) is given by 

-I -1 



D{q,E) = q' 



1 

ap 



rp 



h{q,E) + h{q,E) 



3/2 



h{q,E) 



r(2 + r) 
(l + r)2 



2(1 +r) 
2 + r 



m2E 



(9a) 
(9b) 



where r = mi/ 1712 is the mass ratio. The kernel R"^^ {p,q, E) depends on the quantum 
numbers J^. For positive parity J+ with J > 0, the bound-state amplitude B{p) is a 
two-component column vector and R^^ {p,q, E) is a 2 x 2 matrix with entries 



j+ 



J 1 + r 

+ 



1 + r ■ 2J-1 ^''-^ 
^(J-l)(2J+l)(l + r) 



J 




Q 


Qj-i 






P. 




J-2, 







R 



2J- 1 

= ^J(J + l)(-[l/(l + r) + l + r]Qj 



R: 



j+ 
22 



+{q/p)Qj^i + {p/q)Qj+i), 
J 



1 1 + r 

+ 



l + r ■ 2J + 3 
(J + 2)(2J+l)(l + r; 



(^ + 1) 

Qj+2, 



p^q 
q p^ 



Q 



(lOa) 



(10b) 



(10c) 



2J + 3 

and R21 {p,q, E) = Rf^ {q,P, E). For J > 0, (5j(a;) is a Legendre function of the second 
kind and its argument is x = (1 + r)[p/q + q/p\/2 — m2rE/{pq). For J < 0, Qj = 0. For 
the special case = and for negative parity J~ with J > 0, -B(p) has one component 



and R"^^ is a 1 x 1 matrix. For 0"^, its entry R^{p,q^ E) is given by Eq. (10c). For J , its 



entry is 



R'~ip,q,E) = (1 +r) [Qj+iix) - Qj-iix)] . 



(11) 



These equations are valid for 1/ap < |rp|^/54. 
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IV. RENORMALIZATION GROUP LIMIT CYCLES 



In most channels, the solution to the bound state equation for B{p) decreases rapidly 
enough at large p that the integral in Eq. (|8]) converges as A — )■ oo. In these channels, the 
model is renormalizable: the only parameters that are required are ap and rp. However 
there are also channels in which B{j)) approaches a log-periodic function of p at large p, 
just like in the STM equation. In this case, the integral does not converge, but instead 
approaches a log-periodic function of A as A — oo. If the identical particles are bosons, the 
only channel in which this behavior occurs is 1"*". If the identical particles are fermions, this 
behavior occurs in the channels 0^, l"*", 1~ and 2"*". 

In the field theory framework, the cutoff dependence implies that the model defined by 
the Lagrangian in Eq. ([6]) is not renormalizable in the three-atom sector. In the case of an 
S-wave contact interaction, renormalizability can be restored by adding a three-body contact 
interaction term |30]. In our case with a P-wave interaction, renormalizability can also be 
restored by adding an appropriate three-body contact interaction for each channel in 
which there is log-periodic dependence on A. The renormalization group (RG) flow for its 
coupling constant is governed by an RG limit cycle [Ij. An RG limit cycle is characterized 
by a discrete scaling factor Aq such that every time A increases by a factor of Aq, the coupling 
constant returns to its original value. 

Renormalization can be implemented by eliminating the coupling constant for the three- 
atom interaction in favor of a renormalization scale A^ ^D]- The dependence of any renor- 
malized amplitude on A* can only be log-periodic with discrete scaling factor Aq. The model 
therefore exhibits asymptotic discrete scale invariance at high energies. A discrete scale 
transformation consists of rescaling all momenta and energies by factors of Aq and Aq", re- 
spectively, and also rescaling the parameters as ap X^^^ap and rp ^ -^o^Pj where n is an 
integer. Under a discrete scale transformation, any amplitude changes simply by an overall 
power of Aq. 

The most dramatic consequence of the discrete scale invariance is the Efimov effect. 
Efimov showed in 1970 that if the two-body S-wave scattering length is tuned to infinity, 
there can be an infinite sequence of three-body bound states called Efimov states that have 
an accumulation point at the three-body threshold |3T]. The binding energies of successive 
states differ by the square of a discrete scaling factor Aq that depends on the mass ratios 
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and symmetries. The cases in which the Efimov effect occurs include three identical bosons, 
two identical bosons and a third particle, and two identical fermions and a third particle 
whose mass is smaller than the fermion mass divided by 13.6 (see, e.g. In the case of 
resonant P-wave interactions between one atom and two identical atoms, we have identified 
all the channels in which the Efimov effect occurs. If the identical particles are bosons, 
the Efimov effect occurs only in the 1"^ channel. If the identical particles are fermions, the 
Efimov effect occurs only in the channels 0"^, 1"^, 1~ and 2+. The Efimov effect in the 1~ 
sector was discovered by Macek and Sternberg [221 Our results for the 1~ sector are 
consistent with their results. Unfortunately the Efimov effect in this model is unphysical. In 
the unitary limit, there are negative-probability states at the three-atom threshold, which 
is the accumulation point for the infinitely many Efimov states. 

The discrete scaling factor for each channel can be expressed as Aq = exp(7r/so), where 
isQ is a pure imaginary solution to the equation 



det 



1 



-1 



\J+i 



TT 2J+ 1 



'1 







(12) 



r{2 + r)^ 

and M"^^ is essentially the Mellin-transform of R^^ in Eq. (|8]). For positive parity J+ with 
J > 0, M'^^(s) is a 2 X 2 matrix with entries 

J 1 +r " 



fj-j [i- + 1+] /, 



J-2, 



M- 



J+ 

12 



j+ 

22 



1 + r 2J-1 
(J-l)(2J+l)(l + r 
2J- 1 

v/J(J+l)(-[l/(l + r) + l + r]/j 
+i^fj-i + i-fj+i), 

fj-iJ + 1) [i-^ 



J -I 



J+ 1 1+r 
1 + r ^ 2J + 3 
(J + 2)(2J+l)(l + r) 



J+i 



2J + 3 



fj+2, 



and (s) = M{'2^(— s). The functions fj{s) are 



fAs) 



P.J 



1 + r 



:i3a) 



(13b) 



(13c) 



(14) 



where Pj is a Legendre polynomial if J > and if J < 0, is an operator that shifts the 
argument of fj{s) by ±1, and fo{s) is 

TT sin(s arcsin[l/ (1 + r)]) 



fo{s) 



s cos(s7r/2) 



(15) 
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bosons: 1+ 
fermions: 0+ 
fermions: 1+ 
fermions: 1- 
fermions: 2+ 



r = nij/m^ 



FIG. 1: (Color online) Discrete scaling factor Aq = exp(7r/so) as a function of the mass ratio 
r = mi/m2- If the identical particles are bosons, the Efimov effect occurs only in the = l"*" 
sector. If the identical particles are fermions, the Efimov effect occurs in the 0^, 1^, and 2"*" 
sectors. 

In the special case = O"*" and for negative parity J~ with J > 0, M"^^(s) is a 1 x 1 matrix. 
For 0^, its entry (s) is given by Eq. (13c). For J~ , its entry is 



(16) 



If the identical atoms are bosons, Eq. (12) has a pure imaginary solution isq only for 
= 1+. For equal masses (r = 1), the discrete scaling factor Aq = exp(7r/so) is 116.7. 
If the identical atoms are fermions, Eq. (12) has a pure imaginary solution for O"*", 1"*", 1~, 



J 



and 2+. For equal masses, the discrete scaling factors are 34.2, 239.3, 111.2, and 501.9, 
respectively. The dependence of Aq on the mass ratio r = 1711/7712 for each of the five cases 
is shown in Fig. [TJ The discrete scaling factors all approach 1 as r — )■ 0. 

We proceed to discuss the spectrum of the three-body bound states at rp = predicted 

I I —1/3 

by Eq. (8). In the unitary limit defined by Op = rp = 0, discrete scale invariance 
implies the existence of infinitely many Efimov trimers. They have an accumulation point 
at the three-atom threshold, and the binding energies of successive trimers differ by a factor 
of Aq. As |ap|~^^'^ increases, each of the Efimov trimers disappears through the three- 
atom threshold at a negative ap and through the atom-dimer threshold at a positive ap. 
Discrete scale invariance requires that if an Efimov trimer disappears through a threshold 
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(1/aJ 




FIG. 2: (Color online) Efimov plot showing the energy variable X as a function of with 
rp = for the case of identical fermions, equal masses, and = 0"*". The scaling factor Aq ~ 34.2 
was divided by 17 to fit more states in the plot. 



identical particles 


JP 


So 


Ao 






bosons 


1+ 


0.660 


116.7 


-0.134(2) 


0.452(2) 


fermions 


0+ 


0.889 


34.2 


-0.522(2) 


0.244(2) 


fermions 


1+ 


0.574 


239.3 


-0.25(2) 


0.34(2) 


fermions 


1" 


0.667 


111.2 


-0.576(2) 


0.212(2) 


fermions 


2+ 


0.505 


501.9 


-0.590(1) 


0.19(3) 



TABLE I: The imaginary solutions isq of Eq. ( 12 ) and the discrete scaling factors Aq for the case 



of equal masses mi = m2- The threshold parameters Op and Op for rp = are also given. 

at {ap,rp) = {a*p,r*p), the next deeper Efimov trimer must disappear through the same 
threshold at {X^^a},, Xqv},). 

The bound state spectrum at rp = predicted by Eq. (|8| is illustrated in Fig. |2] for 
the case with identical fermions, equal masses, and quantum numbers O"*". The energy 
variable K = sign{E){fi\E\/ti^y^'^ for several Efimov trimers is shown as a function of ap^^^ 
for rp = 0. An Efimov trimer whose binding energy at unitarity is h'^nl/m disappears 
through the three-atom threshold at Op = — and through the atom-dimer threshold 
at Op = +0.24/k^. The behavior in the other channels is similar. The discrete scaling factors 
for equal masses and the corresponding threshold parameters for all channels are given 
in Table HI 
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V. CONCLUSION 



Whether there is a universal binding mechanism for shallow P-wave states remains an in- 
triguing and important question. Apart from ultracold atoms, such systems occur frequently 
in halo nuclei [3l]. The application of effective field theory to such a system is called halo 
EFT^. 

In this paper, we have discussed the renormalization of the minimal zero-range model for 
P-wave interactions defined by Eq. ([6]). This model is renormalizable in the two-body sector, 
but it has unphysical negative-norm molecules with momentum scale rp. It can be used as 
an effective theory for typical momenta well below the scale rp where the effects of negative 
probability states are suppressed. In the three-body sector, renormalization in some parity 
and angular- momentum channels involves an ultraviolet limit cycle, indicating asymptotic 
discrete scale invariance for large momenta. In the unitary limit a^^^^rp — )■ 0, there is an 
accumulation of three-body states at threshold analog to the Efimov effect. However, this 
limit is unphysical since the negative norm states are driven to threshold and lead to the 
breakdown of the model. 

A promising strategy for dealing with P-wave resonant interactions is to treat the unitarity 
term ik^ from Eq. ([2]) as a perturbation [35j . Such an expansion can be justified rigorously for 
small momenta, k <^ |rp|, and leads to a simplified pole structure of the dimer propagator. 
As an example, we show in Fig. [sjthe dimensionless binding momentum ^JS|Ji\E\/\rp\ as a 
function of the dimensionless cutoff 2A/|rp| for fermions in the 0"'"-channel, with 1/ap = 0, 
rp < 0, and mi/m2 = 1/8. There can be 0, 1, or possibly even 2 three-body bound states for 
dimensionless cutoffs below 1 where the expansion of the dimer propagator can be justified. 
Increasing the mass ratio mi/m2 also increases the critical value of the dimensionless cutoff 
at which the first three-body bound state appears. For equal masses, this critical value is 
larger than 1, so there are no such states without introducing a three-body force. 

The strategy of treating the unitarity term as a perturbation was recently applied to 
^He by Ji, Elster, and Phillips [36]. This halo nucleus consists of an alpha particle and 
two weakly-bound neutrons and can be described as an effective three-body system. The 
alpha-neutron interaction has a strong P-wave resonance while the two neutrons interact 
resonantly in an S-wave channel. Proper renormalization of the halo EFT was achieved by 
introducing an appropriate three-body force that was tuned to give the measured binding 
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FIG. 3: (Color online) The dimensionless binding momentum y^8/i|i?|/|rp| as a function of the 
dimensionless cutoff 2A/|rp| for fermions with = 0^, 1/ap = 0, rp < 0, and mi/m2 = 1/8. 

energy of ^He |36] . Rotureau and van Kolck [371 EH] described ^He in the Gamow shell model 
using a different strategy. They applied a halo effective theory with separable neutron-alpha 
and neutron-neutron interactions and a neutron-neutron-alpha three-body force. The range 
of the separable interactions then acts as an ultraviolet cutoff. The unitarity term can 
be treated non-perturbatively without introducing any negative-probability states. Similar 
strategies have been used for the pionless EFT in the three-nucleon system. 

The successful description of the ^He nucleus within halo EFT [36] suggests that the field 
theory presented in this work can describe three-body phenomena in experiments with ultra- 
cold atoms, provided that the unitarity term in the P-wave scattering amplitude is treated as 
a perturbation. Future applications might include the calculation of the scattering-volume 
dependence of the three-body recombination rate and other few-body reaction rates. 
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